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Abstract 

Background: A few lineages of Group A streptococci (GAS) have been associated with a reemergence of severe 
invasive streptococcal disease in developed countries. However, the majority of the comparisons between invasive 
and non-invasive GAS isolates have been performed for collections of reduced genetic diversity or relied on limited 
typing information to distinguish clones. We characterized by several typing methods and compared a collection of 
160 isolates recovered from normally sterile sites with 320 isolates associated with pharyngitis and recovered in the 
same time period in Portugal. 

Results: Although most of the isolates belonged to clones that were equally prevalent in invasive infections and 
pharyngitis, we identified markers of invasiveness, namely the emm types 1 and 64, and the presence of the speA 
and speJ genes. In contrast, emm4, emmlS, and the ssa and speUM genes were significantly associated with 
pharyngitis. There was a strong agreement between the emm type, the superantigen (SAg) genes and the clusters 
defined by pulsed-field gel electrophoresis (PFGE) profiling. Therefore, combinations of particular emm types and 
SAg genes frequently co-occurred in the same PFGE cluster, but there was no synergistic or antagonistic interaction 
between them in determining invasiveness. Only macrolide-susceptible PFGE clones were significantly associated 
with invasive infections or pharyngitis, while the clones of resistant isolates sharing all other molecular properties 
analyzed were equally prevalent in the two groups of isolates. 

Conclusions: This study confirmed the importance of the widely disseminated emm1-T1-ST28 clone in invasive 
infections but also identified other clones linked to either invasive infections (emm64-ST1 64) or pharyngitis 
(emm4-T4-ST39), which may be more limited in their temporal and geographical spread. Clonal properties like 
some emm types or SAg genes were associated with disease presentation, highlighting the importance of bacterial 
genetic factors to the outcome of GAS infections, although other, yet unidentified factors may also play an 
important role. 
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Background 

Streptococcus pyogenes (Lancefield group A Streptococcus, 
GAS) remains one of the most common human patho- 
gens, being responsible for uncomplicated superficial 
infections of the respiratory tract and skin, such as 
tonsillo-pharyngitis and impetigo, but also causing se- 
vere and rapidly progressing invasive disease such as 
necrotizing fasciitis, bacteremia, streptococcal toxic 
shock syndrome (STSS), puerperal sepsis, pneumonia, 
and meningitis [1]. Although the incidence and severity 
of GAS infections in industrialized countries decreased 
for most of the 20 th century, a reemerge of GAS invasive 
disease has been noted since the late 1980s, both in 
North America and in Europe [2]. The annual incidence of 
GAS invasive disease has been estimated at 2.45/100 000 
for developed countries, with a median case fatality rate of 
15% [3]. 

The increase in the incidence of GAS invasive infec- 
tions has been frequently associated with specific clones, 
raising the possibility that the rise of particularly virulent 
clones was responsible for this reemergence, in particu- 
lar the M1T1 clone which is dominant among invasive 
GAS isolates in most developed countries [4,5]. How- 
ever, a higher representation of a particular clone in in- 
vasive infections may be simply due to a high prevalence 
of that same clone in the general GAS population. To 
address this question several studies have performed 
comparisons between the characteristics of the invasive 
clones and the 5. pyogenes isolates associated with car- 
riage or uncomplicated infections in the same time 
period and geographic region. Some of those studies 
reported a differential distribution of M/emm types or of 
T types between invasive and non-invasive isolates and 
confirmed an association between serotype Tl and Ml/ 
emml and invasive infection [6-9], but many others 
found that the major clones responsible for invasive 
infections had a similar prevalence among non-invasive 
infections [10-12]. However, most of the studies per- 
forming such comparisons were either restricted to 
small numbers of isolates or were limited in the typing 
methodologies used, relying essentially on M/emm 
typing. 

Serotyping of GAS based on protein M, a major sur- 
face virulence factor, has long been used as the gold 
standard for the epidemiological surveillance of the 
infections caused by this pathogen. In recent years it has 
been widely replaced by an equivalent approach based 
on sequencing the hypervariable region of the emm gene 
encoding the M protein. However, recent studies show 
that emm typing alone is not sufficient to unambigu- 
ously identify GAS clones and that it must be comple- 
mented with other typing methods such as pulsed-field 
gel electrophoresis (PFGE) macrorestriction profiling or 
multilocus sequence typing (MLST) [13]. Streptococcal 



superantigens (SAgs) secreted by S. pyogenes play an im- 
portant role in the pathogenesis of the infections caused 
by this species [14]. The profiling of the eleven SAg 
genes described so far (speA, speC, speG, speH, spel, spej, 
speK, speL, speM, ssa, smeZ) can be used as a typing 
methodology [15]. Some studies suggested an association 
between the presence of certain SAg genes or of certain 
SAg gene profiles and invasive infections [10,16], al- 
though others failed to establish such an association, 
reporting instead a strong link between the SAg profile 
and the emm type, regardless of the isolation site 
[12,15]. 

We have previously characterized a collection of 160 
invasive GAS isolates collected throughout Portugal be- 
tween 2000 and 2005, and found a very high genetic di- 
versity among this collection, but with a dominant clone 
representing more than 20% of the isolates, which was 
characterized as emwl-Tl-ST28 and carried the gene 
speA [17]. The aim of the present study was to evaluate 
if the clone distribution among the invasive GAS isolates 
in Portugal reflected the clonal structure of the isolates 
causing pharyngitis, in terms of molecular properties 
and antimicrobial resistance. In order to do that, 320 
non-duplicate isolates collected from pharyngeal exu- 
dates associated with tonsillo-pharyngitis in the same 
time period were studied by emm typing, T typing, SAg 
profiling, PFGE macrorestriction profiling, and selected 
isolates were also submitted to MLST analysis. All iso- 
lates were also tested for their susceptibility to clinically 
and epidemiologically relevant antimicrobial agents. The 
great majority of the clones were found with a similar 
frequency among invasive infections and pharyngitis. 
Still, some clones were shown to have a higher invasive 
disease potential and it was also possible to establish sig- 
nificant associations between some emm types and SAg 
genes and disease presentation. 

Results 

Antimicrobial resistance 

All isolates were fully susceptible to penicillin, quinupristin/ 
dalfopristin, chloramphenicol, vancomycin, linezolid, and 
levofloxacin (Table 1). Among the invasive isolates, 19 (12%) 
were resistant to erythromycin, while the isolates asso- 
ciated with pharyngitis presented significantly higher 
macrolide resistance - 21% (P = 0.016, two-tailed Fisher's 
exact test). Among the invasive macrolide-resistant isolates, 10 
(53%) presented the M phenotype and were therefore 
susceptible to clindamycin, whereas the remaining nine 
(47%) were also constitutively resistant to clindamycin 
(cMLSb phenotype). The proportion of the two pheno- 
types was similar among the pharyngitis isolates, with 
37 isolates (55%) presenting the M phenotype and 30 
(45%) presenting the MLSb phenotype (one with indu- 
cible resistance and the others with constitutive 
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Table 1 PFGE clusters presenting antimicrobial resistant 
isolates collected from invasive infections and 
pharyngitis in Portugal 



PFGE cluster" 


Antimicrobial resistance 6 


No. of resistant isolates 






Invasive Pharyngitis 


c 3S 


Tet 


1 


D 36 


MLS B 


1 




M 


1 


G27 


M 


6 19 




MTet 


1 


H ;?5 


MLS B ,Bac 


6 17 




Tet 


1 


1 


MLS B ,Tet 


1 


J,6 


Tet 


12 1 


K14 


M 


1 


L,3 


MLS B ,Tet 


1 6 




Tet 


2 


M„ 


MLS B ,Tet 


1 


N,o 


Tet 


1 1 




MLS B ,Tet 


1 


0 9 


M 


4 5 


R 6 


M 


3 


s 6 


M 


1 



" Clusters are designated by capital letters and a subscript number indicating 
the number of isolates in each cluster; 

b The antibiotics tested were penicillin quinupristin/dalfopristin, 
chloramphenicol, vancomycin, linezolid, levofloxacin, erythromycin, 
clindamycin, tetracycline, and bacitracin. M, presenting the M phenotype of 
macrolide resistance; MLS B , presenting the MLS B phenotype of macrolide 
resistance; Tet, non-susceptibility to tetracycline; M,Tet, presenting the M 
phenotype of macrolide resistance and resistance to tetracycline; MLS B ,Tet, 
presenting the MLS B phenotype of macrolide resistance and resistance to 
tetracycline; MLS B ,Bac, presenting the MLS B phenotype of macrolide resistance 
and resistance to bacitracin. 

resistance to clindamycin). All the isolates presenting 
the M phenotype of macrolide resistance carried only 
the meJ[A) variant of the mef determinant. The cMLSb 
isolates carried only the ewi(B) gene, except for one 
pharyngitis isolate which also harbored we/(A), and the 
only iMLSs isolate in the collection that presented the 
erm(A) gene. 

In contrast to erythromycin, tetracycline resistance was 
much lower among the pharyngitis isolates when com- 
pared with the invasive group (6% vs 17%, P < 0.001). 
One invasive isolate presented intermediate resistance to 
tetracycline (MIC = 3ug/ml). All the resistant strains car- 
ried the tet(M) gene, except one pharyngitis isolate for 
which no PCR product was obtained for any of the 
screened tetracycline-resistance genes. The tet(V) gene 
was detected in only one pharyngitis isolate, which also 
harbored tet(M), while the genes tet(K) and tet(Qi) were 
not amplified in any of the studied isolates. Overall there 
was a positive association between the genes tet(M) and 
erm{B) (P < 0.001), but that association was not observed 



among the subset of invasive isolates, since only three of 
the 27 tei(M) -positive invasive isolates also carried the 
erm(B) gene (P = 0.178). 

Bacitracin resistance was detected in a total of 23 iso- 
lates (5%), with no significant differences among the two 
types of infection considered. All these isolates expressed 
the cMLSb phenotype of macrolide resistance and were 
tetracycline-susceptible. 

Characterization of GAS clones 

Globally, among the 480 isolates there were 36 emm 
types, 17 T types, and 49 SAg profiles (the genes 
included in each SAg profile are presented in Additional 
file 1). In the subset of 170 isolates (100 from pharyngitis 
and 70 from invasive infections) selected for MLST ana- 
lysis, 49 different STs were identified. Nineteen PFGE 
clusters (groups of > 5 isolates presenting > 80% similar- 
ity on the PFGE profile) were obtained including 268 
pharyngitis isolates and 143 invasive isolates (86% of all 
isolates) (Table 2 and Table 3). Except for R 6 , isolates 
grouped into PFGE clusters presented some variability 
in their emm type, ST, T type, or SAg profile, with most 
variability found in the later two properties. Still, in most 
PFGE clusters the majority of the isolates were charac- 
terized by a single profile of dominant properties. The 
emm diversity among the PFGE clusters differed signifi- 
cantly (Table 4). Within each PFGE cluster, different 
emm types were associated with distinct SAg profiles 
(Table 2 and Table 3), although globally the emm and 
PFGE had a similar predictive power over the SAg pro- 
file (data not shown). 

Unrelated STs within the same PFGE clusters were 
associated with isolates of different emm types, while 
isolates of the same emm type presented the same ST or 
single-locus variants (SLVs) (Table 2 and Table 3). The 
only exceptions were ST39 and ST561 that were both 
associated with cluster G27 and emm% but were double- 
locus variants (DLVs) of each other. In clone I24, four 
distinct STs were found. While ST25 and ST554 were 
SLVs and were both associated with ewm44/61, ST 150 
belonged to a different clonal complex, but was also 
associated with a different emm type (emm75). Finally, 
ST555 despite being associated with an isolate of a dif- 
ferent emm type (emm89) is a SLV of ST25, which may 
explain why this isolate was clustered in I 2 4 and not in 
the major PFGE cluster associated with this emm type 
(C 38 ). 

Isolates expressing the M phenotype of macrolide re- 
sistance belonged mostly to clusters G27, O9, and R6, 
while the majority of MLSb isolates were clustered in 
H 26 and L 13 (Table 1). Isolates belonging to emm types 
1, 4, and 28 were separated into different PFGE clusters 
according to macrolide resistance (B 49 and O9 for 
emml; F 29 and G 2 7 for emw4; H 26 and M n for emm28, 
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Table 2 Properties of the PFGE clusters with >15 GAS isolates collected from invasive infections and tonsillo-pharyngitis 



in Portugal 



PFGE 




No. of isolates (% of total) 
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17 fit 
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3 fit 
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7 (1) 
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6 (3.8) 


16 (5.0) 


2/27/44 (19), NT (2), 12 (1) 
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0 


1 (0.3) 


25(1) 
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0 
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J,6 


64 
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53 


2 (1.3) 


0 


NT (2) 


26 (2) 


11 (D 




74 


0 


1 (0.3) 
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11 (D 


120 (1) 




87 


0 


1 (0.3) 
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38(1) 


62(1) 




89 


0 


1 (0.3) 


B3264 (1) 


43(1) 
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" Clusters are designated by capital letters and a subscript number indicating the number of isolates in each cluster; b NT, non-typeable; c nd, not determined. 



respectively), while the remaining characteristics (T type, 
ST, and SAg profile) of each pair of PFGE clusters were 
the same. 

Bacitracin-resistant isolates were all clustered in PFGE 
H 2 6 and were characterized as ewm28-T28, except for 
one isolate that was emw22-T12. However, this cluster 
was not restricted to bacitracin-resistant isolates, since it 
also included three bacitracin susceptible isolates, two of 
which were also ewm28-T28, while the other was 
emmlS but T non-typable. 

Surface antigen differences between invasive and 
pharyngitis isolates 

The invasive isolates were significantly less diverse than 
the pharyngitis isolates by T typing and SAg profiling 
(Table 5). However, while the emm type distribution 
varied between the invasive and pharyngitis isolates 
(P < 0.001) no differences were noted in the T types. Six- 
teen emm types occurred only in invasive infection or pha- 
ryngitis, but in most cases the small number of isolates 
associated with these emm types prevented the differ- 
ences from reaching statistical significance (Figure 1). In 
contrast, the overrepresentation of emm types 1, 4, 64, 
and 75 in one of the groups was statistically supported. 



SAg differences between invasive and pharyngitis isolates 

A detailed analysis of the SAg profiling results of the iso- 
lates is performed elsewhere [18]. Briefly, the chromoso- 
mally encoded SAg genes smeZ and speG were the most 
frequent among the 480 isolates (n = 461 and 417, re- 
spectively), followed by speC (n = 247), ssa (n = 170), 
spej (n = 157), speA (n = 154), speK (n = 118), speH 
(n = 82), spel (n = 73), and speL and speM, which were 
always detected together (n = 44). The association of 
individual SAg genes with disease presentation was 
tested. In the analysis of these results, the SAg genes 
speG and smeZ were not considered because they were 
present in nearly all isolates, and the genes speL and 
speM were considered as a single entity, since they 
always co-occurred. Individually, the genes speA and spej 
were both associated with invasive isolates (P < 0.001). 
As expected, strains carrying both these genes were also 
associated with invasive infections (P < 0.001), but no 
synergistic effect between the two genes was observed, 
since the presence of one did not significantly increase 
the representation of the other among invasive isolates. 
In contrast, speC (P = 0.002), ssa (P < 0.001), and speL/M 
(P < 0.001) were individually associated with pharyn- 
gitis. The combinations speC+speL/M and ssa+speL/M 
were both associated with pharyngitis (P = 0.004 and 
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Table 3 Properties of the PFGE clusters with <15 GAS isolates collected from invasive infections and tonsillo-pharyngitis 
in Portugal 



PFGE 
cluster" 


emm 
type 


No. of isolates (% of total) 
Invasive Pharyngitis 


T type b 
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22 


0 


1 (0.3) 


12(1) 


21 (1) 


46(1) 


0 9 


1 


4 (2.5) 


5(1.6) 


1 (8), 13 (1) 


10 (9) 


28(4) 


P 8 


78 


4 (2.5) 


4(1.3) 


11 (7), 3/13 (1) 


29 (8) 


409 (3) 


Q 8 


43 


4 (2.5) 


0 


3/13 (2), NT (2) 


11 (4) 


3 (2) 




58 


2 (1.3) 


2 (0.6) 


NT (4) 


17(3), 14(1) 


410 (3), 176 (1) 


R 6 


75 


0 


6(1.9) 


25 (6) 


39 (6) 


150 (2) 


s 6 


9 


1 (0.6) 


4(1.3) 


9(5) 


40(5) 


75 (2) 




12 


0 


1 (0.3) 


12(1) 


33 (1) 


36(1) 



a Clusters are designated by capital letters and a subscript number indicating the number of isolates in each cluster; b NT, non-typeable. 



0.012, respectively), but there was also no synergistic 
effect relative to the presence of a single gene. How- 
ever, the association of speC with pharyngitis isolates 
can be explained by a high frequency of co-occurrence 
of this gene with ssa, since the isolates harboring speC 
without ssa were not significantly associated with any 
of the groups. An interesting situation occurred when 
analyzing the interaction between spej (associated with 
invasive infections) and ssa (associated with pharyn- 
gitis). Among isolates carrying spej, the group that also 
carried ssa was no longer associated with invasive 
infections, while the association of isolates carrying ssa 
with pharyngitis was not significantly altered by the 



Table 4 Simpson's index of diversity and 95% Confidence 
intervals (CI95%) of emm types for each PFGE cluster 



PFGE cluster" 


No. emm types 


SID [Cl 95% ] 


B 49 


2 


0.041 [0-0.118] 


c 38 


2 


0.053 [0-0.151] 


D 36 


2 


0.056 [0-0.159] 


H26 


3 


0.151 [0-0.336] 


'24 


J 


0.163 [0-0.361] 


J,6 


5 


0.533 [0.255-0.812] 


L,3 


5 


0.628 [0.353-0.903] 


N,o 


2 


0.200 [0-0.504] 


Q 8 


2 


0.571 [0.571-0.571] 


s 6 


2 


0.333 [0-0.739] 


" PFGE clusters A 51 , E 30 , 


F 29 , G 27 , K 14 , M,,, 0 9 , P 8 , and R 6 


include only one emm 



type (SID=0). 



presence of spej. This argues for a dominant effect of 
the presence of ssa over that of spej in determining 
the invasive capacity of individual isolates. The associ- 
ation of SAg profiles with disease presentation was 
also tested. Two SAg profiles presented a significant as- 
sociation with invasive isolates, namely SAglO (speA 
+ speG + spersmeZ + ) and SAg46 (speG + smeZ + ) (P < 0.001). 
The remaining profiles were not significantly associated 
with any of the two groups of isolates. 

When the same kind of analysis was performed for 
emm types and individual SAg genes, three combina- 
tions with statistical significance emerged: the associ- 
ation of isolates presenting emml and speA, and emml 
and spej with invasive infections (P < 0.001), and the as- 
sociation of isolates carrying emmTS and speL/M with 
pharyngitis (P = 0.001). In all cases, no synergistic or an- 
tagonistic interaction was detected between emm type 
and SAg gene, since the emm type did not alter the asso- 
ciation of the SAg gene with a particular group of 
isolates. 

Differences between the PFGE clusters found among 
invasive infection and pharyngitis 

The associations described above can be correlated 
with the PFGE clusters which were also different be- 
tween the invasive and pharyngitis groups of isolates 
{P < 0.001), in agreement with the differences found in 
emm types (Figure 1 and Figure 2). All the 19 major PFGE 
clusters occurred in both invasive and pharyngitis iso- 
lates, except for R 6 (eww75-T25-ST150-SAg39), which 
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Table 5 Simpson's index of diversity and 95% Confidence intervals (CI 95 o /o ) of the typing methods used in the analysis 
of the 160 invasive isolates and 320 pharyngitis isolates 



Typing method 


Invasive 




Pharyngitis 






No. types 


SID [Cl 95% ] 


No. types 


SID [Cl 95% ] 


T typing 


13 


0.882 [0.859-0.904] 


17 


0.915 [0.907-0.923] 


emm typing 


30 


0.920 [0.900-0.940] 


26 


0.921 [0.911-0.931] 


PFGE profiling 


30 


0.930 [0.912-0.947] 


44 


0.947 [0.939-0.954] 



SAg profiling 27 0.911 [0.891-0.931] 46 0.941 [0.932-0.951] 



was present only among pharyngeal isolates, but the 
difference did not reach statistical significance due to 
the small number of isolates in this cluster. PFGE dis- 
tinguished several groups of isolates belonging to emm 
types 1 and 4. The difference in the distribution be- 
tween pharyngitis and invasive infection was not found 
for all PFGE clusters containing isolates carrying each 
of these two emm types, but only for those including 
macrolide-susceptible isolates (B 49 associated with inva- 
sive infections and F 2 9 associated with pharyngitis, re- 
spectively, P < 0.001). The emml and emm4 isolates 
expressing macrolide resistance (M phenotype) were 
grouped into PFGE clusters Og and G27, respectively, 
which presented a similar prevalence among invasive 
infections and pharyngitis (Figure 2). PFGE Ji 6 , which 
included all emw64 isolates, was also associated with 
invasive infections (P < 0.001). The emm75 association 
with pharyngitis was not translated into an association 
of a specific PFGE cluster, since the 19 emmTo strains 
were scattered into various PFGE clusters (Table 2 and 
Table 3). 

Not surprisingly, three eww-PFGE cluster combina- 
tions showed significant associations with infection type: 
ewwl-B 4 9 and eww64-Ji 6 were associated with invasive 
infections, while eww4-F 2 9 was associated with pharyngitis 



(P < 0.001). It was not possible to detect any synergistic or 
antagonistic interaction between PFGE and emm type in 
modulating the association of the isolates with either 
group. The same was true for the statistically significant 
combinations between PFGE clusters and individual SAg 
genes, namely the combination of B49 with speA and with 
spej (both associated with invasive infections, P < 0.001), 
and the combination of F 2 g with speC and with ssa (both 
associated with pharyngitis, P < 0.001). 

Discussion 

Several studies yielding conflicting results have attempted 
to compare the clonal composition of GAS populations 
causing invasive and non-invasive infections in order to 
identify particularly virulent clones or properties that may 
be used as epidemiological markers of invasiveness 
[7,8,11,12,16]. However, many of those studies were li- 
mited in the size and diversity of the GAS collections 
studied or in the typing methodologies used, with most of 
them relying essentially on emm typing, which has been 
shown to be insufficient for the complete identification of 
GAS clones [13]. In this work, we used several different 
typing methods to compare a collection of genetically 
diverse GAS isolates recovered from normally sterile sites 
during a period of six years in Portugal [17] with isolates 



P Invasive 
□ Pharyngitis 



T- T- OJ 



emm type 

Figure 1 Distribution of emm types among 160 invasive isolates and 320 pharyngitis isolates. Other includes emm types identified in < 5 
isolates - emm18 (n=4), 25 (n=1), 29 (n=2), 30 (n=1), 43 (n=4), 48 (n=1), 53 (n=3), 58 (n=5), 74 (n=2), 77 (n=4), 90 (n=1), 94 (n=3), 102 (n=3), 103 
(n=1), 113 (n=3), 114 (n=1), 118 (n=1), st106M (n=1), sfG1750 (n=1), srILI 03 (n=1). The asterisk indicates significant differences (P<0.001). 



Friaes et al. BMC Microbiology 201 2, 1 2:280 
http://www.biomedcentral.com/1471 -21 80/1 2/280 



Page 7 of 12 



O 10 



□ Invasive 

□ Pharyngitis 



ffl 



w 



A5I B4g C 38 Dgg E30 F29 G27 H26 I24 Jl6 K14 L13 Mn N10 Og Pg Qg Pig Sg 

PFGE cluster 

Figure 2 PFGE clusters found among 160 invasive isolates and 320 pharyngitis isolates. Approximately 1 1% of invasive and 16% of non- 
invasive isolates were included in PFGE clusters of < 5 isolates that are not represented. The asterisk indicates significant differences (P<0.001). 



recovered from pharyngeal exudates of patients presenting 
with tonsillo-pharyngitis, during the same time period and 
in the same geographical region. The nasopharyngeal mu- 
cosa has been suggested to be the main reservoir for GAS 
isolates associated with invasive infections [19,20]. 

The differences among the GAS clones associated with 
invasive infections and pharyngitis were reflected in anti- 
microbial resistance, with the invasive group of isolates 
presenting lower macrolide resistance and higher tetra- 
cycline resistance, when compared with the pharyngitis 
group. Among isolates belonging to the same emm 
types, namely emml and emw4, only the macrolide- 
susceptible clones were associated with either invasive 
infections or pharyngitis. The macrolide-resistant clones 
of these emm types are reflected in invasive infections 
according to their prevalence in pharyngitis, suggesting 
that these are translating more the antibiotic selective 
pressure than the invasive capacity of the clones. Tetra- 
cycline is not currently used in the treatment of GAS 
infections but resistance to this antimicrobial in S. pyo- 
genes is usually acquired by horizontal transfer, since the 
resistance genes are frequently encoded in mobile gen- 
etic elements with a wide host range [21]. These ele- 
ments often carry macrolide resistance genes as well, 
and in S. pyogenes a significant association between the 
presence of the genes erm(B) and tet{M) has been 
reported and it has been suggested that tetracycline use 
could contribute to the selection of macrolide-resistant 
GAS isolates [21,22]. In our study, the association be- 
tween the presence of the genes erm(B) and tet(M) was 
observed globally, but not among the invasive isolates, 
suggesting that the genetic elements carrying tetracyc- 
line resistance conferring genes may be different be- 
tween the two bacterial populations. 

Bacitracin susceptibility is routinely used for the pre- 
sumptive identification of GAS, although resistant clones 
have been identified in several countries [23-25]. In our 



GAS collection, all the bacitracin-resistant isolates (5%), 
regardless of the type of infection, were clustered in the 
same PFGE clone (H 2 6) and belonged to ST52, although 
one was emw22-T12 while the others were all emm28-T28. 
Isolates with such characteristics had been previously 
reported in Portugal associated with tonsillo-pharyngitis, 
skin infections and asymptomatic carriage [26,27], but 
not with invasive infections. Bacitracin resistance among 
invasive isolates has been previously reported only 
among isolates recovered in France and in San Francisco 
[24,25]. 

Although 74% of the invasive isolates in our collection 
belonged to clones which were equally frequent among 
pharyngitis, suggesting that a significant part of the 
invasive GAS population mirrors the clonal structure 
of the circulating GAS isolates, the remaining isolates 
represented clones that had an enhanced capacity to 
cause invasive disease. We also found significant associa- 
tions between individual properties or pairwise combina- 
tions of properties and disease presentation. Since in 
most cases these were also characteristics of the more 
invasive clones, we cannot exclude that the associations 
of individual properties or pairwise combinations of 
properties can reflect, at least partially, the distribution 
of genetic lineages in the two GAS populations analyzed. 

Individually, emm types 1 and 64 were associated with 
invasive infections. Isolates belonging to these emm 
types presented the only two SAg profiles significantly 
associated with invasive infections. Surprisingly, only 
one of these SAg profiles includes a phage-encoded SAg 
gene {speA). In agreement with our observation, a previ- 
ous study found that within the same PFGE-eww group, 
the SAg profiles significantly associated with invasive 
infections had a smaller number of SAg genes than the 
dominant profiles in pharyngitis [28]. These results sug- 
gest that although some SAg genes may significantly 
contribute to the virulence of S. pyogenes, the rise and 
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success of highly virulent GAS clones may not hinge 
upon the acquisition of phage-encoded SAg genes. Still, 
in our study, the SAg genes speA and spej were both sig- 
nificantly more prevalent among invasive isolates. This 
association was not substantially affected by emm type, 
PFGE clone, nor by the presence of other SAg genes, 
suggesting that speA and spej can be regarded by them- 
selves as markers for invasiveness. Although such associ- 
ation has not been previously noted for spej, the speA gene 
has been frequently associated with invasive infections 
[6,8,16] and the production of SpeA by GAS isolates has 
been linked to streptococcal toxic shock syndrome [29]. 

On the other hand, we identified an association of 
pharyngitis isolates with emm types 4 and 75, and with 
the SAg genes speC, ssa, and speL/M. The association of 
speC with non-invasive infections has been previously 
reported [6,16,30], but in our collection this association 
could be explained simply by a high frequency of co- 
occurrence of this gene with ssa which was strongly asso- 
ciated with pharyngitis, as was also noted in a recent study 
[16]. The presence of the genes speL and speM was not 
previously associated with non-invasive infections. 

Since there is a strong correlation of the SAg profile 
with emm type and of both these properties with PFGE 
type, some of these individual factors frequently co- 
occurred in the same clones. Therefore, combinations of 
these characteristics were also significantly associated 
with disease presentation. However, we could not detect 
any synergistic or antagonistic interactions between most 
of these characteristics, meaning that their co-occurrence 
in a particular isolate does not make it more invasive than 
isolates sharing only one of these characteristics. 

Two PFGE clusters were significantly more prevalent 
among isolates associated with invasive disease than 
among those causing tonsillo-pharyngitis. One of these 
was a cluster of macrolide-susceptible isolates character- 
ized as ewral-Tl-ST28 and by the presence of the SAg 
genes speA, speG, spej, and smeZ (B49), which accounted 
for 18% of the invasive isolates. M1T1 isolates have been 
frequently associated with severe invasive GAS disease, 
and the acquisition of prophage-encoded virulence 
genes, as well as horizontal gene transfer events by hom- 
ologous recombination were implicated in the increased 
virulence of these isolates [31,32]. However, many epi- 
demiological studies found a similar prevalence of this 
clone among invasive and non-invasive isolates [10-12], 
questioning its enhanced invasive capacity. In contrast 
to these findings, but similarly to those of others [6-9], 
we found an association between this clone and invasive 
GAS disease in Portugal, although it can also frequently 
cause milder infections such as pharyngitis (it accounted 
for 6% of the pharyngitis isolates analyzed in this study). 

The other cluster significantly associated with invasive 
infections in Portugal was ] 16 , which was dominated by 



isolates belonging to eraw64-ST164 and carrying the 
SAg genes speG and smeZ. A clone with these character- 
istics has not been previously associated with invasive 
disease and ewm64 has been infrequently reported 
among invasive GAS isolates [4,33,34]. The higher inva- 
sive capacity of this clone cannot be attributed to its 
SAg repertoire, since these isolates do not harbor any of 
the SAg genes associated with invasive infection. Other, 
still unidentified, characteristics may be responsible for 
the properties of this clone. 

In contrast to these PFGE clones, the ¥ 2 g clone of 
macrolide-susceptible isolates characterized by eww4- 
T4-ST39 and harboring the genes speC, ssa and smeZ 
was associated with pharyngitis, suggesting that this 
clone may have a reduced ability to cause invasive dis- 
ease, in agreement with the negative association between 
emm<i and invasive infection that has been suggested 
elsewhere [16]. The association of emmlS with pharyn- 
gitis has not been previously reported and was not trans- 
lated into particular PFGE clusters due to the high 
diversity of emmTo isolates. 

Our data confirms that the widely dispersed M1T1 
clone has enhanced invasiveness but we also identified 
clones with increased or decreased invasive capacity that 
may have emerged locally and that have a more limited 
temporal and geographical spread. The emm alleles and 
the SAg genes characteristic of these clones were asso- 
ciated with particular disease presentations. Other indi- 
vidual emm alleles and SAg genes were also associated 
with a higher propensity to cause invasive infections or 
pharyngitis indicating the importance of these character- 
istics in determining an isolates invasive capacity. 

Other factors that were not evaluated in this study 
may contribute to a different distribution of GAS clones 
in less severe and more severe infections. These include 
bacterial factors, such as the occurrence of mutations in 
transcriptional regulators controlling the expression of 
virulence factors, which seems to play an important role 
in the pathogenesis of some GAS isolates [35]. For other 
clones, the ability to cause invasive infections may be 
more dependent on exploiting host factors, like the HLA 
class II haplotype [36], which may vary in frequency in 
different human populations. 

Conclusions 

This study established links between particular genetic 
lineages and the type of infection, indicating that genetic 
characteristics of the bacteria play an important role in 
determining the outcome of their interaction with the 
human host. The different distribution of clones in the 
two types of infection supports the relevance of PFGE as 
a typing methodology for GAS [13]. This was further 
evidenced by the fact that the macrolide-resistant emm\ 
and emmA PFGE clones were not associated with any 
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particular disease presentation, contrary to the suscep- 
tible clones carrying the same emm types that were 
associated with invasive infections and pharyngitis, res- 
pectively. Moreover, in contrast to other reports 
[12,15] we found associations between particular emm 
alleles and SAg genes and disease presentation. In this 
study, we identified eww4, emmTo, ssa and speL/M as 
independent markers for pharyngitis and emml, 
ewm64, speA, and spej as independent markers for 
invasiveness. Our data re-enforces the multi-factorial 
nature of GAS invasive capacity and highlighted 
lineages and characteristics, in addition to the well 
known M1T1 lineage, that are associated with particu- 
lar disease presentations and that may further increase 
in importance. 

Methods 

Bacterial isolates 

The invasive isolates (w = 160) were collected from nor- 
mally sterile sites, and their partial characterization was 
previously reported [17]. A total of 320 non-duplicate 
GAS isolates were randomly selected among a collec- 
tion of 1604 isolates recovered from pharyngeal exu- 
dates of patients presenting with tonsillo-pharyngitis in 
32 laboratories distributed throughout Portugal, be- 
tween 2000 and 2005, in the proportion of 1:2 (inva- 
sive:pharyngitis) for each studied year. These isolates 
were recovered from pediatric patients (<18 yrs) and 
showed a balanced distribution by gender. The subset 
of macrolide-resistant pharyngeal isolates had been par- 
tially characterized [27,37]. Strains were identified by the 
submitting laboratories and confirmed in our labora- 
tory by colony morphology, p-hemolysis and the pres- 
ence of the characteristic group antigen (Slidex Strepto A, 
BioMerieux, Marcy l'Etoile, France). 

Antimicrobial susceptibility testing 

Susceptibility tests were performed by disk diffusion on 
Mueller-Hinton agar supplemented with 5% defibrinated 
sheep blood, according to the guidelines of the Clinical 
and Laboratory Standards Institute (CLSI) using the fol- 
lowing antibiotic disks (Oxoid, Basingstoke, UK): penicil- 
lin, vancomycin, erythromycin, tetracycline, levofloxacin, 
chloramphenicol, clindamycin, quinupristin/dalfopristin, 
and linezolid. Whenever isolates with intermediate sus- 
ceptibility were identified, the results were confirmed by 
MIC determination using E-test strips (BioMerieux, 
Marcy l'Etoile, France). The macrolide resistance pheno- 
type was determined as previously described [38]. Sus- 
ceptibility to bacitracin was determined for all isolates 
using disks containing 0.05 U of bacitracin (Oxoid, 
Basingstoke, UK), as described elsewhere [27]. 



Detection of antimicrobial resistance genetic 
determinants 

Bacterial DNA was prepared according to the protocols 
of the Centers for Disease Control and Prevention 
(CDC) [39]. Determination of the macrolide resistance 
genotype was performed for strains presenting either the 
M or the MLSb macrolide resistance phenotype, by a 
multiplex PCR reaction with primers to detect the erm 
(B), erm(A) and mef genes, as previously described [40]. 
Isolates carrying the mef gene were subjected to a sec- 
ond PCR reaction in order to discriminate between mef 
(A) and mef[E) [37]. Tetracycline resistant isolates were 
PCR-screened for the presence of the genes tet(K), tet 
(L), tet(M), and tet(0) as previously described [41]. 
Strains harboring each of the resistance genes were used 
as positive controls for the PCR reactions. 

T-typing 

Strains were cultured in Todd-Hewitt broth (Oxoid, 
Basingstoke, UK) at 30°C overnight and treated with 
swine pancreatic extract, using the Auxiliary Reagents 
for Hemolytic Streptococcus Typing (Denka Seiken, 
Tokyo, Japan), and following the manufacturer's instruc- 
tions. T serotypes were determined by slide agglutination 
with 5 polyvalent and 19 monovalent sera (Hemolytic 
Streptococcus Group-A Typing Sera, Denka Seiken). 

emm-typing and SAg gene profiling 

The emm-typing of all isolates was performed according 
to the protocols and recommendations of the CDC, and 
the first 240 bases of each sequence were searched 
against the emm CDC database [39]. Identity of > 95% 
with previously described sequences over the 150 bases 
considered allowed the assignment of an emm type. The 
presence of the SAg genes speA, speC, speG, speH, spel, 
spej, speK, speL, speM, smeZ, and ssa, and of the chro- 
mosomally encoded exotoxin genes speB and speF (used 
as positive control fragments) was assessed in all 160 in- 
vasive and 320 non-invasive GAS isolates by two multi- 
plex PCR reactions as described elsewhere [18]. 

PFGE macrorestriction profiling and MLST 

Agarose plugs of bacterial DNA were prepared as previ- 
ously described [27]. After digestion with Smal or Cfr9I 
(Fermentas, Vilnius, Lithuania), the fragments were 
resolved by PFGE [27]. The isoschizomer Cfr9I was used 
only for the isolates with the M phenotype, which were 
not digested by Smal [13,27]. The macrorestriction pat- 
terns generated were compared using the Bionumerics 
software (Applied Maths, Sint-Martens-Latem, Belgium) 
to create UPGMA (unweighted pair group method with 
arithmetic mean) dendrograms. The Dice similarity coef- 
ficient was used, with optimization and position toler- 
ance settings of 1.0 and 1.5, respectively. PFGE clones 
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were defined as groups of >5 isolates presenting profiles 
with > 80% relatedness on the dendrogram [13]. MLST 
analysis was performed as described elsewhere [42] for 
representatives of each PFGE cluster (a total of 100 non- 
invasive and 70 invasive isolates). When more than one 
emm or T-type was present in the same PFGE cluster, 
isolates expressing different surface antigens were 
selected. Allele and sequence type (ST) identification 
was performed using the S. pyogenes MLST database 
[43]. Whenever new alleles were identified, chromato- 
grams of both strands were submitted to the database 
curator for approval and an allele number was assigned. 
Clonal complexes were determined using the goeBURST 
algorithm implemented in PHYLOViZ [44]. 

Statistical analysis 

The diversities of the different PFGE clusters were com- 
pared using the Simpson's index of diversity (SID) with 
corresponding 95% confidence intervals (Ci95%) [13]. 
Differences in antibiotic resistance between the invasive 
and non-invasive groups of isolates were evaluated using 
Fisher's exact test. P values < 0.05 were considered to in- 
dicate statistical significance. SAg genes, emm types and 
PFGE types were screened for associations with the inva- 
sive group by computing an odds-ratio and an associated 
Fisher's exact test. Additionally, pairs of individual SAg 
genes with each other or with emm types or PFGE types 
were similarly tested for the association of each pairs' 
co-occurrence with the invasive group of isolates. For 
the pairs where at least one of the types individually or 
their co-occurrence were associated (either positively or 
negatively) with the invasive group, two more tests were 
done, to investigate if the association of one of the types 
individually was modified by the co-occurrence of the 
other type in the pair (synergism or antagonism). Con- 
sidering a pair of types A and B, this test compares the 
proportion of invasive isolates among the ones that have 
A type but not B with the same proportion among iso- 
lates that have both A and B types. If the proportions 
are statistically different, according to a Fisher's exact 
test, we can conclude that type B modifies the associ- 
ation of type A with the invasive group of isolates. Con- 
versely, if the proportion of invasive isolates among the 
ones that have the B type but not A differs from the 
same proportion among isolates that have both A and B 
types, type A modifies the association of type B with the 
invasive group. If the isolates that are simultaneously of 
the A and B type show a significantly stronger associ- 
ation with invasive infection than the one observed for 
isolates having either the A or B type, the types are said 
to be synergistic. If, on the other hand, isolates that are 
simultaneously of the A and B type show a significantly 
weaker association with invasive infection than the one 
observed for isolates having either the A or B type, the 



types are said to be antagonistic. All the p-values 
obtained in each step of the screening procedure were 
corrected for multiple testing through the False Discov- 
ery Rate (FDR) linear procedure [45]. 

Additional file 



Additional file 1: SAg genes profiles identified in GAS isolates 
in Portugal. 
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